Quantum Hall Effect on the Hofstadter Butterfly 
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Motivated by recent experimental attempts to detect the Hofstadter butterfly, we numerically 
calculate the Hall conductivity in a modulated two-dimensional electron system with disorder in 
the quantum Hall regime. We identify the critical energies where the states are extended for each 
of butterfly subbands, and obtain the trajectory as a function of the disorder. Remarkably, we 
find that, when the modulation becomes anisotropic, the critical energy branches accompanying a 
change of the Hall conductivity. 



(N 



c3 



-a 

o 
o 



> 

in 

o 
in 
o 
-i— > 

c3 



i 

C 

O 

o 



X 



The problem of a Bloch electron in a magnetic field 
in two dimensions has long been investigated. Theoreti- 
cally, the interplay of the modulation and the magnetic 
field gives rise to a recursive set of energy gaps as a sensi- 
tive function of the field amplitude, which is known to be 
the Hofstadter butterfly Q. It is also predicted that the 
system exhibits the quantum Hall effect when the Fermi 
energy is in each single gap, where the gap-rich structure 
leads to the nontrivial sequence of the quantized Hall 
conductivities Q . On the experimental side, many efforts 
have been devoted to the challenge of observing the Hofs- 
tadter butterfly mainly in lateral superlattices patterned 
on GaAs / AlGaAs heterostructures 0, 0, 13 ■ A remnant 
of the nonmonotonic behavior of the Hall conductivity 
peculiar to the Hofstadter butterfly was observed 

In a two-dimensional (2D) electron system without pe- 
riodic potentials, it is generally believed that the weak 
disorder makes almost all the states localized leaving 
the extended states at the center of the Landau band, 
and the Hall conductivity jumps from one integer to 
another when the Fermi energy sweeps through the ex- 
tended states. It is natural to expect in a periodic system 
that each butterfly subband exhibiting a non-zero Hall 
conductivity has extended states at specific energies. It 
is an intriguing problem how these energies move and 
how the nonmonotonic behavior of the Hall conductiv- 
ity changes as a function of disorder strength, until the 
butterfly structure is eventually destroyed. 

The problems of the quantum Hall effect and Anderson 
localization in the disordered 2D Bloch system have been 
studied by several authors. The evolution of the extended 
states has been investigated in tight-binding lattices, 
from the viewpoint of the disappearance of the quantum 
Hall effect in the weak-field limit 0, S & 03 E3 The 
pair annihilation of extended states with the increase in 
randomness was demonstrated @ and the Hall conduc- 
tivity was also calculated 9]. The nonmonotonic behav- 
ior of the Hall conductivity was demonstrated for the 
Hofstadter butterfly in the presence of disorder [l^. A 
finite-size scaling analysis was performed for a 2D system 
modulated by a weak periodic potential and the critical 
exponent was estimated at the center of the Landau level 
[l3| . A qualitative discussion on the evolution of the ex- 
tended states in the Hofstadter butterfly as a function of 
the disorder for several flux states was given 0, ^| . 

In this paper, we study the scaling property of the 



Hall conductivity a xy in a weakly modulated 2D elec- 
tron system, and identify the critical energies for each 
of butterfly subbands. We obtain the trajectory in the 
energy-disorder space to show how the energies move as 
a function of disorder. We also find that the number 
of critical energies sensitively depends on the anisotropy 
of the periodic potential, which is controllable in a real 
experiment. 

Let us consider a two-dimensional system in a strong 
magnetic field with a square periodic potential 

V p (x,y) = 14 cos— x + VyCos— y. (1) 

The band structure is characterized by the parameter 
<f> = Ba 2 /(h/e), a number of magnetic flux quanta pene- 
trating a unit cell [lj . The disorder potential Vd is taken 
as the randomly distributed delta potential ±uo with 
the concentration 27rZ = 4, where I is the magnetic 
length and n.; is the number of scatterers in a unit area 
[lfij . The energy scale for the disorder is then given by 
The magnetic field is assumed to 
be sufficiently strong so that V p and T are smaller than 
the cyclotron energy hw c , and so that we can consider 
only the lowest Landau level. To deal with a finite sys- 
tem, we consider the larger unit cell L x L with L = Ala 
(M: integer) and apply the periodic boundary condition. 
We calculate the Hall conductivity using the Kubo for- 
mula for zero temperature, 



iL 2 



(2) 

where e a is the energy of the eigenstate |a), Ef the Fermi 
energy, and v t the velocity operator. We take into ac- 
count the mixing of the states between Landau levels 
n = and 1 up to the first order of V/(hui c ) with V 
being V p or V d . 

We first consider an isotropic modulation, V x — V y . 
The calculations are performed for two different fluxes, 
= 3 and 3/2. In both cases, the lowest Landau level 
splits into three subbands in the absence of the disorder, 
while the distributions of the Hall conductivity a xy are 
different as shown in Fig. ^ Specifically, a xy for three 
subbands becomes (0,1,0) (in units of —e 2 /h) in <j) = 3 
and (1, —LI) in <t> — 3/2, which satisfies the Diophantine 
equation 0. 
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FIG. 1: Plots of the Hall conductivity a xy (in units of — e 2 /h) 
for a nondisordered electron in the modulation V x /V y = 1 
with 4> — 3 (upper) and 3/2 (lower). Dashed lines represent 
the density of states in units of l/(V y a 2 ). 
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FIG. 2: (a) Hall conductivity a xy (in units of —e 2 /h) cal- 
culated for disordered systems with V x /V y = 1, F/Vy = 0.5, 
= 3, and L/a = 8, 12, and 16. The dashed line represents 
the density of states in units of l/(V y a 2 ). (b) Relative values 



of a xy measured from the smallest [L/a 



sample, (c) 



Inverse localization length (in units of a) estimated from the 
Thouless number. The vertical dashed line penetrating the 
panels represents the critical energy. 



Figures [21 and |2| show the numerical results calculated 
for the disordered systems with several sizes for = 3 
and 3/2, respectively. The panel (a) shows o~ xy with the 
density of states, and (b) the difference in o~ xy measured 
from the smallest sample. We also show in the panel 
(c) the inverse localization length l/L\ oc estimated by 
the Thouless number method [Hi- We divide the energy 
axis into 200 columns and every quantity is averaged over 
a number of different samples for each column. 
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FIG. 3: Plots similar to Fig. for V x /V y = 1, V /V v = 0.5, 
and = 3/2 with L/a = 12, 16, 20, and 24. The panel (d) 
shows the trajectories of the critical energies as a function of 
r. Numbers represent the corresponding Hall conductivities. 



We can clearly see that the subbands with zero Hall 
conductivity (the upper and lower subbands in (f> — 3) are 
all localized as naturally expected, and the Hall conduc- 
tivity there rapidly reaches the quantized value. Each of 
other bands contains a critical energy where the localiza- 
tion length diverges (l/Li oc = 0). The Hall conductivity 
for these delocalized bands has fixed points at a xy = 1/2, 
which agree quite well with the critical energies estimated 
from Li oc . a xy off the fixed points always approaches 1 in 
the area a xy > 1/2 and in the area <r xy < 1/2, leading 
to the Hall plateau in an infinite system (shown as step 
like lines in (a)). Therefore, it is natural to define the po- 
sitions of the extended states as the points of a xy = 1/2, 
which is exactly the method adopted in the following. 
The panel (d) in Fig. [3] shows the trace of the critical 
energies identified as the positions of a xy — 1/2 for sev- 
eral values of r. Three lines with the Hall conductivities 
1,-1, and 1 become closer as the magnitude of disorder 
increases, and combine into one critical energy with the 
Hall conductivity 1. 

The localization property is strongly affected by the 
configuration of the periodic potential itself. Here, we 
investigate how the above results are altered when we 
make the potential anisotropic (i.e., V x ^ V y ). Figure 01 
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FIG. 4: (a) Hall conductivity a xy (in units of — e /h) cal- 
culated for disordered anisotropic systems with V x /V y = 0.2, 
Y/Vy = 0.25, <f> = 3/2 and L/a = 16, 24, 32, 40. The dashed 
line represents the density of states in units of l/(V y a 2 ). (b) 
Relative values of a xy measured from the smallest (L/a = 16) 
sample, (c) Trajectories of the critical energies as a function 
of r. Numbers represent the corresponding Hall conductivi- 
ties 



FIG. 6: (a) Traces of the critical energies for several values of 
Vx/V y . For V x /V y = 0.4, we present the energy regions where 
the error bar of the Hall conductivity reaches a xy = 1/2. 
(b) Phase diagram showing the number of critical energies in 
the disorder (r/V r ll )-anisotropy (V x /V y ) space. '?' indicates 
indeterminate. The dashed line shows 4Ea/V y with Eg being 
the subband gap in the ideal system without disorder. 
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FIG. 5: Plots of the Hall conductivity a xy (in units of 
—e 2 /h) for a nondisordered electron in anisotropic systems 
with Vx/Vy = 0.5, 0.2, 0.05 and and (j> = 3/2. The dashed 
line represents the density of states in units of l/(V y a 2 ). 



shows the Hall conductivity calculated for the disordered 
system with V x /V y = 0.2 and <j> = 3/2. Remarkably, 
we found five fixed points at a xy = 1/2, and they likely 
lead to the critical energies with the Hall conductivities 
1, — 1, 1, — 1, and 1 in an infinite system, shown as a zig- 
zag line in (a). We also calculated the Thouless number, 
but failed to estimate the localization length around the 
extended region, since we cannot resolve the exponent of 
the size dependence of the Thouless number. However, 
the robust fixed points and the monotonic behavior with 
L in the Hall conductivity strongly suggests the existence 
of the critical energies at those points. The trajectories 
of the fixed points for several values of T are shown in 
Fig. ^c). At a certain point, the upper and lower pairs 
annihilate and the center is left, which is topologically 
different from the isotropic case. 

The origin of the five critical energies can be explained 
as follows. Figure 13 shows the Hall conductivity a xy (E) 
for clean systems with the fixed flux = 3/2 and vari- 
ous ratios of V x /V y . The gaps between subbands rapidly 
shrink as V x /V y goes apart from 1 and disappear at 
V x = (ID modulation), where a xy (E) behaves just 
monotonically from to 1. When a very small V x is 
present, we see that a xy (E) is deviated from that in 
V x = only around the subband gaps (E — E gap £V X ). 
There the perturbative V x mixes the states in V x = to 
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produce positive and negative Hall currents so that the 
Hall conductivity is quantized in the gap. If a small but 
finite disorder is introduced into this system, the criti- 
cal energies with opposite Hall conductivities should ex- 
ist correspondingly on each side of the gap, as long as 
r is smaller than the gap width and the Hall plateau 
survives. The electron-hole symmetry requires another 
critical energy right at the Landau band center. When 
V x is switched off, the critical energy pairs at the gaps 
annihilate and the central one is left. 

We show the corresponding plots for various values of 
V x /V y in Fig. Efa), to explain how the isotropic (three 
critical energies) and anisotropic cases (five energies) are 
connected. We see that the center branch in the isotropic 
case splits into three branches at V x /V y « 0.4. The points 
at r = show the energies of a xy = 1/2 in clean systems. 
The critical energies should converge to these points in 
the limit r — ► 0, as long as a xy = 1/2 is the critical point 
for the disordered systems with a finite size. The number 
of such points changes from 3 to 5 at V x /V y = (V x /V y )o w 
0.36, which defines the turning point in the evolution of 
the critical energy. The results (although not included in 
Fig. EJa)) show that, for V x /V y = 0.35 < (V x /V y ) , there 
are five critical energies when T is sufficiently small, and 
for Vx/Vy — 0.45, there are only three critical energies. 
In the case V x /V y — 0.4, however, it is difficult to resolve 
the critical energies around E = 0, partly because of the 
numerical accuracy of the calculated a xy . Therefore, we 
can present only the energy regions where the error bar 
of the Hall conductivity reaches a xy = 1/2. 

We give the phase diagram in Fig. Elb), showing 
the number of critical energies in the disorder (T/V y )- 



anisotropy (V x /V y ) space. The newly found five-energy 
phase occupies a considerable area in the region of high 
anisotropy and low disorder. The dashed line shows 
AEo/V y as a function of V x /V y , where Eg is the subband 
gap in the ideal system without disorder. The critical dis- 
order where the number of critical energies changes from 
five to one or three to one roughly follows the curve as 
might be expected. 

While we concentrated on the flux = 3/2 here, it is 
expected that, in other fluxes, the butterfly gaps simi- 
larly shrink as the modulation goes to ID and the pair 
annihilation of the extended states occurs at the butter- 
fly gaps. We also assume that such a behavior does not 
basically change in higher Landau levels. We expect that 
controlling the QHE by the anisotropy proposed here is 
sufficiently promising for a real experiment, if a modu- 
lation is given by the external gate electrodes as in the 
recent experiment . The existence of the fixed point at 
&x y = 1/2 is compatible with the the two-parameter scal- 
ing theory 1 1 8| . and investigating the (a xy , a xx ) scaling 
in the periodic system is an important future work. 
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